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ABSTRACT 

A  complex  proportionate-type  normalized  least  mean  square  algo¬ 
rithm  is  derived  by  minimizing  the  second  norm  of  the  weighted 
difference  between  the  current  estimate  of  the  impulse  response  and 
the  estimate  at  the  next  time  step  under  the  constraint  that  the  adap¬ 
tive  filter  a  posteriori  output  is  equal  to  the  measured  output.  The 
weighting  function  is  assumed  positive  but  otherwise  arbitrary  and 
it  is  directly  related  to  the  update  gains.  No  assumptions  regarding 
the  input  signal  are  made  during  the  derivation.  Different  weights 
(i.e.,  gains)  are  used  for  real  and  imaginary  parts  of  the  estimated 
impulse  response.  After  additional  assumptions  special  cases  of  the 
algorithm  are  obtained:  the  algorithm  with  one  gain  per  impulse  re¬ 
sponse  coefficient  and  the  algorithm  with  lower  computational  com¬ 
plexity.  The  learning  curves  of  the  algorithms  are  compared  for  sev¬ 
eral  standard  gain  assignment  laws  for  white  and  colored  input.  It 
was  demonstrated  that,  in  general,  the  algorithms  with  separate  gains 
for  real  and  imaginary  parts  have  faster  convergence. 

Index  Terms —  Adaptive  filtering,  convergence,  least  mean 
square  algorithms. 

1.  MOTIVATION 

The  complex  least  mean  square  (CLMS)  adaptive  filter  [1]  was  orig¬ 
inally  proposed  to  extend  the  least  mean  square  algorithm  from  real¬ 
valued  signals  to  complex-valued  signals.  To  date  there  exists  no 
equivalent  extension  of  proportionate-type  normalized  least  mean 
square  algorithms  (PtNLMS)  from  real-valued  signals  to  complex- 
valued  signals.  PtNLMS  algorithms  are  a  family  of  adaptive  algo¬ 
rithms  which  update  the  estimate  of  the  impulse  response  in  a  man¬ 
ner  which  is  proportional  to  the  magnitude  of  the  current  estimate  of 
the  impulse  response.  Examples  of  PtNLMS  algorithms  include  the 
proportionate  normalized  least  mean  square  (PNLMS)  [2]  algorithm, 
and  /r-PNLMS  (MPNLMS)  [3]  algorithm. 

In  this  paper,  the  PtNLMS  algorithm  is  extended  from  real- valued 
signals  to  complex-valued  signals.  The  resulting  algorithm  is  named 
the  complex  PtNLMS  (CPtNLMS)  algorithm.  The  CPtNLMS  algo¬ 
rithm  updates  the  real  and  imaginary  parts  of  the  unknown  impulse 
response  by  applying  separate  real-valued  gains.  Additionally,  sev¬ 
eral  simplifications  of  the  CPtNLMS  will  be  proposed. 

This  paper  is  organized  in  the  following  manner.  First,  nota¬ 
tion  and  the  adaptive  filter  framework  will  be  introduced.  Next  the 
derivation  of  the  CPtNLMS  algorithm  will  be  presented  as  well  as 
several  simplifications  and  interpretations.  After  that,  examples  of 
CPtNLMS  algorithms  and  their  implementations  will  be  discussed. 
Next  simulation  results  for  various  CPtNLMS  algorithms  will  be 


presented  and  finally  a  conclusion  summarizing  the  results  found  in 
this  paper  will  be  given. 

2.  COMPLEX  ADAPTIVE  FILTER  ALGORITHM 
INTRODUCTION  AND  NOTATION 

2.1.  Notation 

Let  A  denote  an  arbitrary  complex-valued  matrix.  The  conjugate, 
transpose,  and  conjugate  transpose  of  the  matrix  A  are  denoted  by 
A*,  AT,  and  AH\  respectively.  Additionally,  Re(A)  and  Im(A) 
represent  the  real  and  imaginary  parts  of  the  complex  matrix  A. 
These  operators  have  the  same  meanings  when  employed  on  vectors 
and  scalars  where  applicable.  Let  the  ith  component  of  any  vector 
a  be  denoted  as  Oi  and  the  (i,  j)^  entry  of  any  matrix  A  as  [A]jj 
throughout  this  work.  Next,  for  vector  a  with  length  L  we  define  the 
function  Diag(a)  as  an  L  x  L  matrix  whose  diagonal  entries  are 
the  L  elements  of  a  and  all  other  entries  are  zero.  Finally,  we  define 
3  =  \/-T- 

2.2.  Complex  Adaptive  Filter  Framework 

All  signals  are  complex  throughout  this  work.  Let  us  assume  there  is 
some  complex  input  signal  denoted  as  x(fc)  for  time  k  that  excites  an 
unknown  system  with  complex  impulse  response  w.  Let  the  output 
of  the  system  be  y(k)  =  x.H  (k)-w(k)  where  x(fc)  =  [x(k),x(k  — 
1), . . . ,  x(k—L+l)]T  and  L  is  the  length  of  the  filter.  The  measured 
output  of  the  system,  d(k),  contains  complex-valued  measurement 
noise  v(k)  and  is  equal  to  the  sum  of  y(k)  and  v{k).  The  impulse 
response  of  the  system  is  estimated  with  the  adaptive  filter  coefficient 
vector,  w (fc),  which  has  length  L  also.  The  output  of  the  adaptive 
filter  is  given  by  y(k)  =  xw(fc)w(fc).  The  error  signal  e(k)  is 
equal  to  difference  of  the  measured  output,  d(k)  and  the  output  of 
the  adaptive  filter  y(k). 

For  notational  convenience  we  will  suppress  the  time-indexing 
notation  and  instead  use  the  following  convention.  For  an  arbitrary 
vector  a  we  denote  a (fc  +  1)  =  a+  and  a (k)  =  a. 

3.  COMPLEX  PROPORTIONATE-TYPE  NORMALIZED 
LEAST  MEAN  SQUARE  ALGORITHM  DERIVATION 

3.1.  Complex  PtNLMS  Algorithm  Using  Separate  Gains  for  Co¬ 
efficient  Imaginary  and  Real  Parts 

The  estimated  weight  vector  can  be  written  in  terms  of  real  and  imag¬ 
inary  parts  w  =  w_r  +  j'wi,  where  w_r  and  w/  are  vectors  of  length 
L  consisting  of  the  real  and  imaginary  components  of  w,  respec¬ 
tively.  Motivated  by  the  derivation  of  the  normalized  least  mean 


978-1  -467 3 -0046-9/ 1 2/$26 .00  ©2012  IEEE 


3285 


ICASSP  2012 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

MAY  2012 


2.  REPORT  TYPE 

N/A 


3.  DATES  COVERED 


4.  TITLE  AND  SUBTITLE 

Complex  Proportionate-Type  Normalized  Least  Mean  Square  Algorithms 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


6.  AUTHOR(S) 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 
5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Naval  Research  Laboratory  Radar  Division  Washington,  DC  20375,  USA  REPORT  number 

9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADA561051.  AOARD-CSP-111007  International  Conference  on  Acoustics,  Speech  and  Signal 
Processing  (37th)  (ICASSP  2012)  Held  in  Kyoto,  Japan  on  March  25-30,  2012.  U.S.  Government  or  Federal 
Purpose  Rights  License.,  The  original  document  contains  color  images. 

14.  ABSTRACT 

A  complex  proportionate-type  normalized  least  mean  square  algorithm  is  derived  by  minimizing  the 
second  norm  of  the  weighted  difference  between  the  current  estimate  of  the  impulse  response  and  the 
estimate  at  the  next  time  step  under  the  constraint  that  the  adaptive  filter  a  posteriori  output  is  equal  to  the 
measured  output.  The  weighting  function  is  assumed  positive  but  otherwise  arbitrary  and  it  is  directly 
related  to  the  update  gains.  No  assumptions  regarding  the  input  signal  are  made  during  the  derivation. 
Different  weights  (i.e.,  gains)  are  used  for  real  and  imaginary  parts  of  the  estimated  impulse  response. 

After  additional  assumptions  special  cases  of  the  algorithm  are  obtained:  the  algorithm  with  one  gain  per 
impulse  response  coefficient  and  the  algorithm  with  lower  computational  complexity.  The  learning  curves 
of  the  algorithms  are  compared  for  several  standard  gain  assignment  laws  for  white  and  colored  input.  It 
was  demonstrated  that,  in  general,  the  algorithms  with  separate  gains  for  real  and  imaginary  parts  have 
faster  convergence. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

SAR 

4 

19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


square  (NLMS)  algorithm  [4]  let  us  consider  the  following  mini- 


mization  problem 

min 

w+ 

(Wfl  -  Wfl)TMi?1(wJ  -  Wfl) 

+(wj  —  w/)tM71(w+  —  w  i 

such  that 

A  H  -  + 

a  =  x  w 

where  Mb  and  M/  are  real-valued,  nonnegative,  diagonal  matrices, 
which  have  the  property  that  Tr[Mfl]  =  Tr[Mj]  =  L. 

The  method  of  Lagrange  multipliers  will  be  used  to  cast  this 
constrained  minimization  problem  into  one  of  unconstrained  mini¬ 
mization.  Prior  to  performing  this  step  the  following  definitions  are 
introduced  to  aide  in  the  derivation.  Let 


UJ  = 

w  R 

M  = 

W I 

Mfl  0 

0  M/ 


(2) 


The  following  form  of  xHw  will  also  be  employed 


T  •»  |  T  ~  ,  •  T  ~  T 

xrw R  +X,WJ  +  jxHwr  -  jx7  w R 


=  [Xfl,  xf] 


r  H  ■  H  i 

[x  ,  JX  ] 


Wfl 

W I 

w  R 

W  / 


+  j[-X/.  Xfl] 


wH 

W  l 


r  H  ■  H  i  - 
[X  ,  JX  ]w. 


(3) 


Using  these  definitions  the  minimization  problem  can  be  rewritten 
as 


min  J{u>+) 

uy^~ 


=  (u>+  -  w)TM_1(w+  -  u>) 

+  A  (d-  [xH,jxw]a>) 

+  (d*  ~  [xT,  -jxT]d>)  . 


(4) 


Returning  our  attention  to  (6)  this  equation  can  be  rewritten  as 


2(d) +  -  d>) 


Mrx* 

jMjx’ 


+  A* 


Max 

-jM/x 


(9) 


Using  the  definition  in  (2)  allows  us  to  write 

2(wr  -  wr)  =  AMflX*  +  A*Mrx 

2(w^"  —  wj)  =  jAM/x*  —  jA*M/x.  (10) 

Next  using  (10)  and  (8)  we  form 

w+  —  w  =  (wr  —  wr)  +  ji'wf  —  w  /) 

=  8x*  +  yMrrx  -  ^M,x*  +  ^M;x 

=  ^  (M_r  -  Mr)x*  +  ^-(Mr  +  Mr)x 

=  [xg(Ma  +  MJ)xe-lx^(Ma-MJ)xe]  _ 

[xff(Ms  +  M/)x]2  —  |xt(Mr  —  M/)x|2  V 

[xfl(Mfl  +  Mf)xe  -  xfl(Mfl  -  M;)x*e*l 
[x»(Mfl  +  Mr)x]2-|xT(Mfl-Mr)x|2  V  ' 

(ID 


Next  the  step-size  parameter  /3  is  introduced  to  allow  control  over 
the  update.  The  resulting  CPtNLMS  algorithm  with  separate  gains 
for  real  and  imaginary  coefficients  is  given  by 


+P 

+P 


=  w 

[x"(Mr  +  Mj)xe*  -  xt(Mr  -  Mr)xe] 
[xfl(Ms  +  Mr)x]2  -  |xt(Mr  -  M/)x|2 

[xff(Ms  +  M/)xe  -  xs(Mii  -  M,)x*e* 
[x»(Mr  +  M,)x]2  -  |xT(M R  -  Mr)x|2 


(Mn  -  M,)x* 

-  (M_r  +  M/)x. 

(12) 


3.2.  CPtNLMS  Algorithm  Simplifications 


Next  taking  the  derivative  of  J(u>+)  with  respect  to  d>+  and  setting 
the  result  to  zero  yields 


9J(t2>+) 

du}+ 


2M_1(d>+ 

0. 


■  x* 

x 

* 

-A* 

.  ix 

.  _^X  . 

(5) 


The  CPtNLMS  algorithm  given  in  (12)  can  be  simplified  under  cer¬ 
tain  conditions.  For  instance  it  is  straightforward  to  show  that  by 
setting  Mfl  =  Mj  =  G  in  (12)  the  algorithm  is  reduced  to 


w  +  /3 


Gxe 

xHGx 


(13) 


As  it  turns  out  (13)  is  the  solution  to  the  following  minimization 


Multiplying  (5)  from  the  left  by  M  and  rearranging  terms  allows  us 
to  write 


u>) 

=  AM 

x’ 

+  A*M 

X 

.  7X* 

.  -3*  . 

Next  we  employ  the  relationships  given  by  [xH ,  jxH]  (d>+  — d>)  =  e 
and  [xT ,  —  jxT]  (d)+  —  d>)  =  e*  to  form  the  following  linear  system 
of  equations 


* 

e 

xh(Mr  +  M;)x  xt(M_r  —  M/)x 

'  A  ' 

e 

xff(Mfl-Mi)x*  xff(Mfl  +  M7)x  _ 

A* 

It  is  straightforward  to  invert  the  matrix  given  in  (7)  which  allows  us 
to  find 


x  _  2  xg(Mfl  +  Mj)xe*  -  xT(Mfl  -  M/)xe 
[x«(Mfl  +  M,)x]2  -  |xT(Mfl  -  M7)x|2  ' 


(8) 


min 


such  that 


(w+  —  w)wG  1(w+  —  w) 

i  H  *  ~\~ 

a  =  x  w 


(14) 


where  G  is  the  stepsize  control  matrix.  This  version  of  the  CPtNLMS 
algorithm  uses  one  real-valued  gain  to  simultaneously  update  both 
the  real  and  imaginary  parts  of  the  estimated  coefficient. 

A  separate  simplification  to  (12)  can  be  made  if  it  is  assumed 

that 

xT(Mfl  —  M/)x  ~  0.  (15) 

Using  this  approximation  (12)  is  reduced  to 


=  w  + 


Mfl(xe  +  x*e*)  +  M/(xe  —  x*e*) 
x-ff(Mfl  +  Mj)x 


■  (16) 


This  version  of  the  CPtNLMS  algorithm  reduces  the  computational 
complexity  of  ( 1 2)  and  will  be  referred  to  as  the  simplified  CPtNLMS 
algorithm. 
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Table  1 


Stepsize  Control  Matrix  for  CPtNLMS  Algorithm 
Using  One  Real- valued  Gain  per  Coefficient 


F[\wi(k)\,k] 

=  Specified  by  the  user 

T min  (&) 

=  P  max{5p,  F[|u)i(fc)|,  fc], . . . , 

F[\wL(k)\,k]} 

7 i(k) 

=  max{7  min(k),F[\wi(k)\,k]} 

9i(k ) 

_  71 M 

i£f-i7i« 

G(k) 

=  Diag{c/i(fc), . . . ,  <?r,(fc)} 

3.3.  Alternative  Representations  of  the  CPtNLMS  Algorithm 

There  are  alternatives  ways  to  represent  (11).  For  instance  one  could 
rearrange  terms  such  that 

w+  —  w  =  (wj  —  w  R)  +  j{  wf  —  W/) 

„  A*x  +  Ax*  „  A*x  —  Ax* 

=  M  r - - - +  Mr - - - . 

(17) 


Using  (8)  it  is  easy  to  show  that 

Re  [ex]  —  Re  ^ 


A*x  +  Ax* 


A*x  —  Ax* 


=  2- 


B(Mg+M,)x 


ex 


: 


cH(Mfi  +  M/)x- 
Im  [ex]  +  Im  |j 


T(Mjj— Mj)x|2 
=h(Mh+M7)x 

’(Mfl-M,)x 
f(Ma+M,)xU 


^(M«  +  Mt)x-^MK. 


(18) 


Combining  (  17  )  and  (18),  and  introducing  the  stepsize  parameter  /3 
will  result  in  an  alternate  form  of  the  CPtNLMS  algorithm.  In  this 
version  of  the  CPtNLMS  algorithm,  the  estimated  impulse  response 
coefficients  are  updated  by  one  purely  real  term  multiplying  Mn  and 
another  purely  imaginary  term  multiplying  M/. 


4.  CPtNLMS  ALGORITHM  AND  IMPLEMENTATIONS 

In  this  section  several  examples  of  CPtNLMS  algorithms  are  pre¬ 
sented.  The  stepsize  control  matrix  for  the  CPtNLMS  algorithm  us¬ 
ing  one  real-valued  gain  per  coefficient  is  shown  in  Table  I.  The  term 
F[|,u)z(fc)|,  fc],  with  l  £  {1,  2, . . . ,  L},  governs  how  each  coefficient 
is  updated.  In  the  case  when  F[|uij(fc)|,  fc]  is  less  than  7 min,  the 
quantity  7 min  is  used  to  set  the  minimum  gain  a  coefficient  can  re¬ 
ceive.  The  constant  Sp ,  where  fip  >  0,  is  important  in  the  beginning 
of  learning  when  all  of  the  coefficients  are  zero  and  together  with 
p,  where  p  >  0,  prevents  the  very  small  coefficients  from  stalling. 
G  (fc)  =  Diag  {gi(k),  . . . ,  <71,  (fc)}  is  the  time-varying  stepsize  con¬ 
trol  diagonal  matrix.  The  constant  S  is  typically  a  small  positive 
number  used  to  avoid  overflowing.  Some  common  examples  of  the 
term  F[\uii(k)\,k\  are  F[\wi(k)\,  k]  =  1,  F[\wi(k)\,k]  =  \uii(k)\, 
and  F[|'u)z(fc)|,  k]  =  ln(l  +  p\wi(k)\),  which  result  in  the  com¬ 
plex  NLMS  (CNLMS),  complex  PNLMS  (CPNLMS),  and  complex 
MPNLMS  (CMPNLMS)  algorithms,  respectively.  The  parameter  p 
is  used  to  adjust  the  amount  of  compression  created  by  the  logarithm 
function. 

Similarly  the  CPtNLMS  algorithm  using  separate  real-valued 
gains  for  coefficients  real  and  imaginary  parts  is  shown  in  Table  II. 
In  this  case  the  user  specifies  two  functions  FR[\wRii(k)\,k]  and 
F/[|u)/!;(fc)|,  k]  which  govern  how  the  real  and  imaginary  parts  of 
a  coefficient  are  updated.  Next  a  minimum  gain  for  the  real  and 
imaginary  parts  are  calculated,  7 min,R(k)  and  7 min,i(k),  respec¬ 
tively.  Finally,  the  stepsize  control  matrix  for  the  real,  Mh((),  and 


Table  II 


Stepsize  Control  Matrix  for  CPtNLMS  Algorithm 
Using  Separate  Real- valued  Gains  for 
Coefficient  Real  and  Imaginary  Parts 


FR[\wRti(k)\,k] 

=  Specified  by  the  user 

F/[ \wpi(k)\,  k] 

=  Specified  by  the  user 

=  pma\{5p,FR[\wRp(k)\,k], . . . , 

FR[\wR,L{k)\,k]} 

r)min}I  (&) 

=  pmax{5p,  Fj[  ru/,i(fc)  ,  fc], . . . , 

Fi[\wip(k)\,  fc]} 

lR,l(k) 

=  max{7 min,R(k),  FR[\wR}i(k)\,k]} 

7/,  i{k) 

=  max{7 mtn,/(fc),F/[  w/,i(fc)  ,fc]} 

mR,i(k) 

7fi,l  W 

z  £it=l  7R,i(fc) 

mi,i(k ) 

71,1  (fc) 

=  Diag{mfl,i(fc), . . .  ,771^,1,^)} 

M  i(k) 

=  Diag{77i/,i(fc), . . .  ,m/,z,(fc)} 

imaginary,  M /(fc),  parts  are  obtained.  An  example  of  the  terms 
FR[\wR,i(k)\,k]  and  Fr[|u)/,i(fc)|,  fc]  is  FR[\wR}i(k)\,  k]  = 
|thfi,i(fc)|  and  Fi[\wi,i(k)\,  k]  —  (fc)| -  This  yields  the  com¬ 

plex  PNLMS  algorithm  using  separate  gains  for  coefficient  real  and 
imaginary  parts.  In  a  similar  vein  the  complex  MPNLMS  algorithm 
can  be  created  by  setting  FR[\wR,i(k)\,  k]  =  ln(l  +  /i|uii{!i(fc)|)  and 
Fi[\wi^i(k)\,  k]  =  ln(l  +  p\wpi(k)\). 

5.  SIMULATION  RESULTS 

In  this  section  we  compare  the  mean  square  error  (MSE)  versus  iter¬ 
ation  for  several  CPtNLMS  algorithms.  In  Figures  1  and  2  we  plot 
the  real  and  imaginary  parts  of  the  impulse  response.  The  impulse 
response  used  has  length  L  =  512.  The  stepsize  control  parame¬ 
ter  /3  was  set  to  a  value  of  0.5  for  all  of  the  CPtNLMS  algorithms. 
The  input  signal  used  in  the  first  set  of  simulations  was  white,  cir¬ 
cular,  complex,  Gaussian,  stationary  with  power  cr%  =  1.  The  noise 
used  in  all  simulations  was  white,  circular,  complex,  Gaussian,  sta¬ 
tionary  with  power  a 1  =  10-4.  The  parameter  values  p  =  0.01, 
5  —  0.0001,  Sp  =  0.01,  and  p  =  2960  were  also  used. 

The  MSE  learning  curve  performance  for  the  CNLMS,  CPN¬ 
LMS  with  one  real- valued  gain  per  coefficient,  simplified  CPNLMS 
with  separate  gains  for  real  and  imaginary  parts,  CPNLMS  with  sep¬ 
arate  gains  for  the  real  and  imaginary  parts,  CMPNLMS  with  one 
real-valued  gain  per  coefficient,  simplified  CMPNLMS  with  separte 
gains  for  the  real  and  imaginary  parts,  and  CMPNLMS  with  sepa¬ 
rate  gains  for  the  real  and  imaginary  parts  algorithms  are  depicted 
in  Figure  3.  Here  we  see  that  the  CPtNLMS  algorithms  outperform 
the  CNLMS  algorithm  in  the  transient  regime.  The  simplified  algo¬ 
rithms  perform  virtually  the  same  as  the  exact  separate  gain  algo¬ 
rithms.  Additionally,  the  separate  gain  for  real  and  imaginary  parts 
of  coefficients  version  result  in  better  convergence  than  the  single 
real-valued  gain  versions. 

In  Figure  4  the  MSE  performance  is  displayed  for  a  colored  in¬ 
put  signal.  The  input  signal  consists  of  colored  noise  generated  by  a 
single  pole  system  as  follows: 

x(k)  =  yx(k  —  1)  +  a(k)  (19) 

where  x(0)  =  a(0),  a(fc)  is  a  white,  circular,  complex,  Gaussian, 
stationary  process  with  power  =  1,  and  complex  pole  7.  The 
value  7  =  0.3560  —  0.8266j  was  used  in  this  simulation.  The  mag¬ 
nitude  of  7  is  equal  to  0.9,  which  implies  a^.  =  <t^/(1  —  1 71 2)  = 
5.263.  Similar  trends  hold  for  the  colored  and  white  input  signal 
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Fig.  1.  Real  Part  of  Impulse  Response 

IMPULSE  RESPONSE 


Fig.  3.  CPtNLMS  Learning  Curve  Comparison 


Ip  ; 

LEARNING  CURVE  COMPARISON 


CNLMS 

CPNLMS  (One  Gain) 
CPNLMS  (Simplified) 
CPNLMS  (Separate  Gains) 
CMPNLMS  (One  Gain) 
CMPNLMS  (Simplified) 
CMPNLMS  (Separate  Gains) 


50  100  150  200  250  300  350  400  450  500 

COEFFICIENT 


2000  4000  6000  8000  10000  12000  14000  16000  18000 
ITERATIONS 


Fig.  2.  Imaginary  Part  of  Impulse  Response 


Fig.  4.  CPtNLMS  Learning  Curve  Comparison  with  Colored  Input 


cases.  Of  course  the  convergence  rate  is  lower  because  of  the  col¬ 
ored  input. 

6.  CONCLUSION 

In  this  work  we  have  introduced  a  family  of  complex  PtNLMS  al¬ 
gorithms  which  were  derived  by  minimizing  the  second  norm  of  the 
weighted  difference  between  the  estimated  impulse  response  at  the 
current  time  step  and  the  next  time  step  under  the  constraint  that 
the  adaptive  filter  output  at  the  next  time  step  equals  the  measured 
output  at  the  current  time.  The  choice  of  the  stepsize  control  ma¬ 
trix  is  arbitrary  during  the  derivation  which  allows  the  user  to  design 
algorithms  to  suit  their  needs.  Additionally,  it  was  shown  through 
simulation  that  adapting  the  real  and  imaginary  components  of  the 
estimated  impulse  response  independently  of  each  other  results  in  in¬ 
creased  MSE  rate  of  convergence.  Moreover,  it  can  be  noticed  from 
the  simulations  that  the  performances  of  the  complex  algorithm  with 


separate  gains  for  real  and  imaginary  parts  and  its  simplified  version 
are  nearly  indistinguishable. 
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